Can low-frequency breathers exist in a quasi-lD crystal? 
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We investigate a quasi- ID crystal: 2D system of coupled linear chains of particles with strong intra-chain 
and weak inter-chain interactions. Nonlinear dynamics of one of these chains when the rest of them are fixed 
is reduced to the well known Frenkel-Kontorova (FK) model. Its continuum limit, sine-Gordon (sG) equa- 
tion, predicts two types of soliton solutions: topological solitons and breathers. It is known that the quasi- ID 
topological solitons exist also in a 2D system of coupled chains and even in a 3D model of a polymer crystal. 
Numerical simulation shows that the breathers inherent to the FK model do not exist in the system of chains. 
The effect changes scenario of kink-antikink collision at small velocities: it always results only in intensive 
phonon radiation while kink-antikink recombination in the FK model results in long-living low-frequency sG 
breather creation. The reason of the difference is that the model of coupled chains catches 'acoustic' part of 
phonon spectrum of a crystal while the FK model does not. Low-frequency SG breathers in a crystal intensively 
emit resonant 'acoustic' phonons and come to ruin. 



The Frenkel-Kontorova (FK) model (a linear chain of har- 
monically coupled particles on the sine substrate) [1] is the 
the most commonly used and comprehensively investigated 
(see monographs [2, 3]) ID model of a crystal. In the case 
of weak substrate potential (weakly discrete system) it seems 
to be especially appropriate for polymer crystals: quasi- ID 
topological soliton-like excitations predicted by its continuum 
limit, sine-Gordon (sG) equation, exist in a 2D system of cou- 
pled chains [4, 5] and even in a 3D model of a polymer crys- 
tal (see, for example, [6]). SG equation is the only nonlinear 
wave equation of type 



u tt 



](u) = 



(1) 



which possesses also one-parametric family of exact solutions 
in the form of low-frequency breathers [7]. Frequencies of the 
breathers fill the gap between zero and the minimal frequency 
in phonon spectrum fi^(O). If the breather frequency ap- 
proaches Qfk(0), the breather amplitude tends to zero, and 
the breather width - to infinity ('phonon' limit). If the breather 
frequency tends to zero, the breather approaches a full kink- 
antikink profile. 

Numerical simulations [8, 9] show that in the FK model 
the sG breathers survive, and, although lose the energy due 
to resonances of odd multiples to the breather frequency with 
phonon frequencies, have their lifetime long enough even in 
the case of strong discreteness (more than hundred periods 
when the third harmonic to its frequency becomes higher than 
the upper phonon band edge [8]). In the case of weak discrete- 
ness losses of energy are hardly perceptible [9] . 

In connection with studying real physical quasi- ID systems 
such as long Josephson junctions and quasi- ID ferromagnets 
there emerged many works treating behavior of sG breathers 
under action of perturbations breaking exact integrability: dis- 
sipative and diverse conservative terms (see [10, 11], and [12] 
and references therein). Analytical treatment of the problem 
is possible as perturbation in the inverse scattering transform 
[10] or as multiple- scale asymptotic expansion [11] in the 
limit of high breather frequencies, and one can obtain some 
estimates in general case [12]. As one can easily predict, the 
breather lifetime proved to be long if perturbation is small. 



In nonintegrable models with substrate potentials sufficiently 
different from the sine function breather- like long-living non- 
linear excitations are observed numerically ((/> 4 - [13], double 
sG, square well potential - [14]). For the </> 4 model it is shown 
[15] that the radiation rate of a small- amplitude 'breather' lies 
beyond all orders in asymptotic expansion. 

All this allows one to look on such breathers as being 'el- 
ementary excitations' in a crystal, together with kinks and 
antikins (topological solitons) and phonons. This implies 
that the breathers can noticeably contribute to thermodynamic 
properties of a crystal [16] and even must be used in phe- 
nomenological approaches to sG thermodynamics instead of 
phonons [17, 18]. We show that this conclusion based on anal- 
ysis of the FK model (a chain on a substrate) is not valid for a 
more realistic model of a crystal: a system of coupled chains. 

Let us take a system of coupled linear chains of (classi- 
cal) particles (fig.l). To catch the main physical meaning of 
the model it is enough to allow inter-chain interactions only 
between particles of the nearest neighboring chains. Then 
Hamiltonian of the system is written as 
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+ U ( r rn,n;j)}, (2) 



where the dot denotes time derivative, Co is the period of a sep- 
arate chain, c - the longitudinal period of the crystal, u m ^ n is 
longitudinal deviation of the particle (m, n) from its equilib- 
rium position (shown in fig. 1) in the crystal (we keep transver- 
sal deviation y m ,n = 0), the potential U (r myn; j) describes in- 
teraction of the n-th particle in the ra-th chain with the (n+j)- 



th particle in the (m - 
between the particles 



l)-th chain, 



m,n;j 



being the distance 
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The ground state of the system (u 
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FIG. 1: 2D model of quasi- ID system: intra-chain interactions are 
much stronger than inter-chain ones (weak discreteness limit). We 
choose unit length so that c = 1 in the equilibrium ground state of 
the crystal. 



energy per particle 

E(b,c,c ) = (c-c ) 2 + U ( R J^ 
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where Rj = [6 2 + c 2 (j + l/2) 2 ] 1 / 2 is the distance between the 
n-th and the (n + j)-th particles of the ra-th (m is odd) and 
(m + l)-th chains. Equilibrium values of b and c minimize the 
expression at given Co . Equivalently, if one chooses c as unit 
length, one can find equilibrium values of cq and b. Hereafter 
we imply that the crystal is initially in its equilibrium ground 
state with c = 1. 

The present 2D model of quasi- ID crystal was first intro- 
duced in [4] . One can take into account also transversal dis- 
placements of particles [5]. The model allows existence and 
propagation of quasi- ID topological soliton-like excitations. 
In [4, 5] have been used the Morse potential of particle inter- 
actions U (r m>n; j) as very suitable for numerical calculations. 
Here we exploit the more physical Lennard- Jones potential 
(truncated): 
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where the truncation function f(r) = {1 — tanh[/i(r — do)]}/ 2 
(/i ~ 1 ? do ^> rg) is introduced for convenience of numerical 
calculations. It allows one to avoid taking into account inter- 
actions of the particles placed one from another farther than 

r « d . 
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FIG. 2: Low-frequency sG breather in the crystal (e = 0.007, 
M = 51, N = 400): oscillations in the FK model (the case of 
immobile neighboring chains) (a) and degradation in the full model 
of interacting mobile chains (b). We show the displacements U26, n 
of the 26th chain containing the breather in successive time moments 
t. Breather frequency is cob — 0.0157. 



It is known [19] that if the equilibrium distance between 
particles ro falls into the interval 0.91 < ro < oo shape of 
the substrate potential is close to the sine function. We have 
chosen ro = 1.67 (do = 20, /i = 2) because it corresponds to 
model of polyethylene crystal with 'united atoms' [6, 20]. At 
this value of ro the substrate generated by immobile neighbors 
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V(u)=2 [U(r J (u))-U(R J )i 
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where rj(u) = [b 2 + (u + j + 1/2) 2 ] 1 / 2 (note that Rj = 
rj(0)), is the sine function accurate within 0.1%: V(u) ~ 
e[l — cos(2ttu)], where substrate amplitude e = 0.1757e. So 
dynamics of one chain when the rest of them are fixed is re- 
duced to the well known FK model: 

H FK = + jK+i - u n ) 2 + V(u n )}. (5) 



The width of a static kink in the model of polyethylene crys- 
tal with 'united atoms' (about 30 periods) coincides with the 
width of a static kink in our model of coupled chains if the 
intensity of inter-chain interactions e = 0.0007. We have also 
considered the cases of stronger interactions e = 0.007, 0.07. 
The first two cases correspond to limit of weak discereteness. 
In the last case the sound velocity in transversal direction is 
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equal to one in longitudinal direction (see table I). When the 
chains are assembled into the crystal the transversal equilib- 
rium period appears to be b = 1.5666 independent on e. Only 
Co is ^-dependent. 

The system of equations of motion for the quasi- ID crystal 
takes the form: 

•• - dH <n 

n = 0,±1,±2, m = 0,±1,±2, ... 

with the Hamiltonian (2). In numerical simulations we con- 
sidered the dynamics of a bounded rectangular fragment of 
the crystal (1 < n < TV, 1 < m < M) with fixed boundary 
conditions in both directions. 

We have compared behavior of a sG breather in the FK 
model (all the chains are kept immobile except one - with 
number m = (M + l)/2 - containing the breather) and in 
the model of coupled chains (all the chains are mobile). 

Numerical simulation shows that the breather in the FK 
model enjoys regular stable oscillations for a very long time 
(fig. 2 (a)). The situation changes drastically if we allow all 
the chains in the crystal to move. The breather quickly comes 
to ruin. Its lifetime is less than two its periods (fig. 2 (b)). The 
effect is observed by all three values of the inter-chain inter- 
actions e = 0.07, 0.007, 0.0007. The destruction results from 
the intensive emission of phonons into the neighboring chains 
(see fig. 3). Energy of the breather spreads to all the particles. 
So one can conclude that the low-frequency sG breathers are 
absent in the model of coupled chains. 

With this in mind, one can suppose the difference between 
the two models under study in scenario of kink-antikink re- 
combination when they collide with small velocities. Indeed, 
in the FK model (with the sine potential as well as with the 
double sine or square well ones) kink and antikink can form 
oscillating breather-like state, their energy remaining for a 
long time localized [14] - see fig. 4 (a), while in the model 
of coupled chains their energy scatters at ones with phonons - 
see fig.4 (b). 

The reason of the effect observed seems to be that the 
two models of a crystal possess qualitatively different phonon 
spectra. Namely, the FK model has a gap between uj = and 
the lower edge of the spectrum ^fk(O) while in the model of 
coupled chains - like in a genuine crystal - the minimal pos- 
sible frequency is uJmin — 0- Let us show it. As we have 
chosen particles numeration not coinciding with one based on 
translation of the crystal cell, phonon modes have the more 



TABLE I: Dependence of the longitudinal s x and transversal s y 
sound velocities, characteristic frequencies uj m ax and ^x(O) on 
the value of the parameter e. 
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Sx Sy/b COmax 


£W(0) 


0.07 


0.7530 0.4916 2.1174 


0.6952 


0.007 


0.9781 0.1555 2.0120 


0.2197 


0.0007 


0.9978 0.0492 2.0012 


0.0694 




FIG. 3: State of the crystal (M = 51, N = 400) of coupled 
chains (e = 0.007) resulting from destruction of low-frequency 
(ujb = 0.0157) sG breather placed onto the 26th chain. The dis- 
placements Um,n are shown in the time moment t = 150. 
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FIG. 4: Kink-antikink recombination in a chain of the crystal (e — 
0.007) having immobile (a) and mobile (b) neighboring chains. Col- 
lision of solitons in the FK model (soliton velocity s = 0.005) results 
in creation of a low-frequency sG breather. Collision of solitons in 
the model of all moving chains (soliton velocity s = 0.25) results in 
intensive radiation of phonons. 

complicated form: 

U2m,n = Aexpi[qin + q22m — u;t], 
U2m+i,n = Aexpi[gi(n - 1/2) + q 2 (2m + 1) - ut],(l) 
n = 0,±1,±2, m = 0,±1,±2, ... . 

where 4 « 1, and gi, q 2 G [0,7r]. Substituting the anzatz 
(7) into the linearized system of equations (6) with imposed 
periodic boundary conditions in both directions one can obtain 
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Now the minimal frequency is Qfk(0) = (^^^Kj) 1 / 2 > 0- 
Dependence of the velocities s x , s y , and the frequencies Umax 
and ^fx(O) on the parameter of inter-chain interactions e is 
presented in the Table I. 



So, with such a phonon spectrum, the model of interact- 
ing chains can not possess even approximate low-frequency 
breather solutions. Their frequencies fall into phonon band, 
the resonance interaction between a breather and a phonon 
takes place and the breather energy is transmitted to the 
phonons, resulting in quick breather degradation. And evi- 
dently it is the case in real (for example, polymer) crystals. 



FIG. 5: Dispersion surface to — Cl(qij qz) for the model of coupled 
chains. The model parameters ro = 1.67, e — 0.07 . The curve 
on the surface is the dispersion curve for corresponding FK model 
(approximation of immobile chains) uo = 7r/2). 

the dispersion equation 

ft(gi,tf2) {2(1 - cos ^1 ) 
+4^i^-[l - cos((j + l/2)gi) cos^]} 1 / 2 , (8) 

3=0 

where rigidities are Kj = U"(Rj)(j + l/2) 2 /R 2 j + 
U' \Rj)b 2 1 'Rj. Values of the rigidities are in direct propor- 
tion to the parameter of inter-chain interation e\ for e = 0.07 
they are K = 0.199, K x = -0.061, K 2 = -0.014, 
Ks = —0.002. We have presented the plot of the disper- 
sion surface for e = 0.07 in fig. 5 together with the dispersion 
curve for the corresponding FK model. 

Dispersion equation (8) gives the minimal u m in = 
£1(0, 0) =0 and the maximal 

+oo 

OOmax = ft(7T, q 2 ) = 2(1 +^K j ) 1 ' 2 

j=o 

possible frequencies, the velocity of longitudinal long 
phonons 

+ CO 

s x = lim n( ? i,0)/gi = [1 + V (j + 1/2)%] 1/2 , 

j = -oo 

and the velocity of transversal long phonons 

+ CO 

s y = b lim n(0,q 2 )/q 2 = b( V Kj) 1 ' 2 . 

q2^0 z — ' 
j=-oo 

The cut of the dispersion surface (8) at q 2 = tt/2 produces the 
dispersion curve for the corresponding FK model: 

+00 

n FK (q) =n(q, tt/2) = {2[1 - cos( 9 )] + 2 Kj} 1 ' 2 . 

j=-oo 
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